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1 INTRODUCTION 



! The integrated Sachs-Wolfe (iSW) effec t JSachs & Wolfd 



2002 



iRees & Sciama 1968'; Hu & Sugi vam3 Il994 ICooraT 
iSchiifer & Bartelinanri 2006), which refers to the frequency 
change of cosmic microwave background (CMB) photons if 
they cross time evolving gravitational potentials, is a direct 
probe of da rk energy because i t vanishes in cosmologies with 
n,„ = 1 dCrittenden & Turokl [l996). By now, it has been 
detected with high significance with a number of different 
tracer objects jposalba et alj I2OO3I: iBoughn & Crittendenl |2004 



Nolta et al.' '2004'; 'Padmanabhan et alJ I2OO5I; 'Oiannantonio et al.' 
2006; Pietrobon et al. 2006; Gaztanaga et al. 2006; Cabre et al. 
2006; Vielva et al. 2006; Rassat et alj|2007l ; iMcEwen et alj|2007l ; 
Giannantonio et al...2008i) . 



ABSTRACT 

The subject of this paper is a quantification of the impact of uncertainties in bias and bias 
evolution on the interpretation of measurements of the integrated Sachs-Wolfe effect, in par- 
ticular on the estimation of cosmological parameters. We carry out a Fisher-matrix analysis 
for quantifying the degeneracies between the parameters of a dark energy cosmology and 
bias evolution, for the combination of the PLANCK microwave sky survey with the EUCLID 
main galaxy sample, where bias evolution b(a) - + {\ — a)ba is modelled with two pa- 
rameters Zjq and ba- Using a realistic bias model introduces a characteristic suppression of the 
iSW-spectrum on large angular scales, due to the altered distance-weighting functions. The 
errors in estimating cosmological parameters if the data with evolving bias is interpreted in the 
framework of cosmologies with constant bias is quantified in an extended Fisher-formalism. 
We find that the best-fit values of all parameters are shifted by an amount comparable to the 
statistical accuracy: The estimation bias in units of the statistical accuracy amounts to 1.19 
for Q.,„, 0.27 for erg, and 0.72 for w for bias evolution with ba - 1. Leaving ba open as a free 
parameter deteriorates the statistical accuracy, in particular on f2,„ and w. 
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the bias undergoes a slow evaluation, which is linked to the lin- 
ear growth function, and decreases slowly towards unity. Numeri- 
cal in v estigation s have found equivalent results (e.g. iBlanton et al] 
I2OOOI) . IHo et alj ( I2OO8I) apply iSW-tomography to measurements 
of the cross-spectrum Cry{{) and the tracer autocorrelation Cyyil) 
using families of tracer objects peaking at different redshifts for 
controlling uncertainties in the tracer redshift distribution, aiming 
at disinguishing flat and curved cosmologies with the iSW-effect. 
Revisiting the iSW-det ection in the NVSS radio source catalogue, 
iRaccanelli et"ai] ( I2OO8I) point out effects of bias evolution on the 
interpreation of the iSW-signal. In their case the observational bias 
(due to luminosity evolution) can correct the discrepancies for the 
iSW-prediction for a ACDM model which arose by using an up- 
dated redshift distribution for the tracer objects. 



1967 



Up to now, the aforementioned studies assumed a constant 
bias in the interpretation of the signal, although it is well estab- 
lished by theoretical studies that the bias of the tracer objects 
evolves by as much as 50% from redshift unity to today. Bias and in 
particular bias evolution has bee n the t opic o f a number of papers, 
both a nalyt ical and numerical. IftvI (Il996h . iTegmark & Peebles! 
(11998') and Basila kos & Phonij bOOlh looked at a two compo- 
nent fluid composed of dark matter and galaxies, and derived mod- 
els for bias evolution in perturbation theory: In their calculations 
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We set out to perform a Fisher-analysis to quantify the magni- 
tude of varying tracer bias on iSW-measurements and to investigate 
the basic degeneracies between the cosmological model, and in par- 
ticular the dark energy equation of state properties, and cosmolog- 
ical bias evolution. We provide a quantification of degeneracies for 
the combination of the PLANCK CMB-observation with the EU- 
CLID galaxy sample. In a further step, we investigate the error in 
parameter estimation if the data with evolving bias is interpreted 
in terms of a cosmology with constant bias, u sing the extension o f 
the Fisher-matrix formalism worked out by jTaburet et al.ll200^ . 
with focus on the estimation of the dark energy equation of state 
parameter w. 
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After compiling the key formulce describing structure forma- 
tion in dark energy cosmologies in Sect. (2] we introduce the two 
relevant observational channels in Sect.[3]and derive the spectrum 
Cry{{) between iSW-temperature perturbation t and the galaxy den- 
sity y in Sect. |4] We quantify the degeneracies between the cos- 
mological parameters using a Fisher-matrix analysis in Sect. [5] and 
extend this formalism to describe the parameter estimation bias re- 
sulting from interpreting the iSW-signal with bias evolution in a 
model with constant bias in Sect. [6] A summary of our results is 
given in Sect.|7] 

As cosmologies, we consider the family of spatially flat ho- 
mogeneous dark energy models with constant dark energy equa- 
tion of state, and with Gaussian adiabatic initial conditions in the 
cold dark matter field. Specific parameter choices for the fiducial 
model in the Fisher-matrix analysis are Hq = lOOAkm/i/Mpc with 
h = 0.72, 0,„ = 0.25, o-g = 0.8, w = -1.0 and = 1, with a 
non-evolving bias parameter of unity. 



2 COSMOLOGY AND STRUCTURE FORMATION 
2.1 Dark energy cosmologies 

In a spatially flat dark energy cosmology with the matter density 
parameter n„, the Hubble function H(a) = dlna/df is given by 



+ (1 - n,„)exp 



' 1 + w{a) 
aa 



(1) 



with the dark energy equation of state w(a). The value w = -1 
corresponds to the cosmological constant A. The relation between 
comoving distance x (given in terms of the Hubble distance dn = 
c/Hq) and scale factor a is given by 



(2) 



with the speed of light c. The dark energy equation of state w(a) can 
be parameterised b y its first order Taylor expansion with respect 
to the scale-factor a l iChevallier & Polarskil200lllLinder & JenkinsI 
l2003h . 



w(a) = Wo + (1 - ajWa- 



(3) 



The conformal time, which is related to the cosmic time t by the 
differential drj = dt/a, follows directly from the definition of the 
Hubble function. 



da 



1 



a^H{a) ' 



(4) 



in units of the Hubble time = 1 /Hq. Hence, rj is defined in com- 
plete analogy to the comoving distance x- 



2.2 CDM power spectrum 

Inflationary models suggest that the CDM power spectrum P(k), 
which describes the fluctuation statistics of the Gaussian density 
field, {S(k)6(k')} = (27r)^(5o(A:+A:')^'(fc) in the case of homogeneous 
and isotropic fluctuations, can be written 



P{k) <x k'-'T^k), 

with the transfer function ( lBardeenet^ll986h . 
ln(l -I- 2.34^) 



(5) 



T(q) 



■ (l + 3.89^ + (\6Aqf + {5A6qf + {6.1\qf) ' 



where the wave vector q is given in units of the shape parameter F 
il^h. P(k) is normalised to the value erg on the scale R = 8 Mpc/h, 



dkk^W^(kR)P{k), 



(6) 



with a Fourier-transformed spherical top-hat W(x) = 3ji(x)/x as 
the filter function. denotes the spherical Besse l function of 
the first kind of order £ jAbramowitz & Stegunlll973) . 



2.3 Structure growth in dark energy cosmologies 

The homogeneous growth of the overdensity field, 6(x,a) = 
D+(a)d(x, 1) is described by the growt h function D+(a), which is 
a solution to the differential equation ( Wang & Steinhardill998l : 
Turner & White 1997; Linder & lenkins 2003), 



d- 1 / dlnH\ d 3 

— 0+ + - 3 + — — = —n„(a)D+(a). 

da-^ a\ dlna / da 2a^ 



(7) 



In the standard cold dark matter (SCDM) cosmology with Q.,„ = 1 
and 3 -I- d In /f/d In a = |, the solution for D+{a) is simply the scale 
factor itself, D+(a) = a. This solution determines the initial con- 
ditions D+(0) = and d/daD+(0) = 1, due to early-time matter 
domination and causes the iSW-effect to vanish in n,„ = 1 cos- 
mologies. 



2.4 Bias and bias evolution 

Dark matter haloes hosting galaxies form by spherical collapse 
at peaks in the cosmic density field. The fractional perturbation 
An/{n) in the spatial number density n of ga laxies is related to th e 
overdensity S = Ap/p by the bias parameter jBardeen et aljl983) . 



An _ ^Ap 
(n) p ■ 



(8) 



The bias parameter b is in general scale dependent and subjected 
to a time-e volut ion, as it slowly decreases w ith time towards unity. 
lFrvl l ll996l) and iTegmark & PeeblesI i 1998h consider the cosmolo- 
giecal evolution of bias in perturbation theory for small perturba- 
tions An/n « 1 and find 



D+(a) - I + ho 



(9) 



with the growth function D+(a) and the value bo of the bias pa- 
rameter today. Motivated by their findings, we use a simple, linear 
model for bias evolution with two parameters bo and b„. 



b(a) = bo + (I - a)ba. 



(10) 



in order to accomodate effects from bias evolu tion if 
the g alaxy density is strong ly perturbed I B Ian ton et al.l 
120001; iBasilakos & PliorfisI 1200 ll) and from velocity bias 
I Percival & Schafeil 20081) . In the course of structure forma- 
tion, thebias of any object tends to unity, bo = 1, as first shown by 
IftvI (11996'). which we will assume applies to the galaxy sample. 
By comparing the observed galaxy correlation function in the 
VMOS-VLT survey with the da rk matter correlation function 
derived from n-body simulations, iMarinoni et al.l ( l2005h propose 
the approximation to the bias evolution law. 



b(z) = 1 + (0.03 ± 0.01)(1 +zy 



(11) 



which if linearised for values of a in the vicinity of 1, yields as a 
lower limit a value of b^ - 0.1. 
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3 LINE OF SIGHT EXPRESSIONS 

3.1 iSW-efTect and large-scale structure tracers 

The iSW-effect is caused by gravitational interaction of a CMB 
photon with a time-evolving potential O. The frac tional perturba- 
tion T of the CM B temperature Tcmb is given by l lSachs & Wolfd 
ll967l:lRees feS ciama 1968) 

T = = - d77 — = -- Axa^Hia) — , (12) 

J CMB t J ori Jo oa 

where rj denotes the conformal time. In the last step, we have re- 
placed the integration variable by the comoving distance x, which 
is related to the conformal time by dx = -cdrj = -cdt/a, and the 
time derivative of the growth function has been rewritten in terms 
of the scale factor a, using the definition of the Hubble function 
da/dt = aH(a), with the cosmic time t. The gravitational potential 
<b follows from the Poisson equation in the comoving frame, 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

scale factor a 



AO = AnGa-p„,{a) 6 



(13) 



where Newton's constant G is replaced with the critical density 
Pcnt(a) = 3H-(a)/(SnG),pm(a) = n„,(a)p^nt(a). 



3HgQ.„, 

2a 



(14) 



Substitution yields a line of sight expression for the linear iSW- 
effect T (integrating along a straight line and using the flat-sky ap- 
proximation), sourced by the linear density field 5, 



3n„, r^" , d /Z)+\ A-' 



(15) 



with the inverse (dimensionless) Laplace operator A ' /d^ solving 
for the (dimensionless) potential if, 

A-' 

V = -ir^, (16) 

with the Hubble distance = c/Hq. Due to the achromaticity of 
the iSW-effect, a measurement of the temperature perturbations of 
the microwave sky can not distinguish between primary and sec- 
ondary, iSW-induced anisotropics. For that reason, the iSW-effect is 
measured in cross-correlation with a tracer of the large-scale struc- 
ture. 



Figure 1. Time evolution of the source terms: /)+(«) for the density field 
(solid line), d(D+/a)/da of the iSW-eflect (dashed line) and Ah(a) = 
hinio) - bo of the linear biasing model (dash-dotted, from top to bottom 
for bo = 1.5, 1.4, 1.3, 1.2, 1.1), with ACDM as the cosmological model. 



with j8 = 3/2 and zo = 0.64, which results in a median redshift of 
Zmed = 0.9. We define the mean bias in a bias evolution model 

biz). 



dz p(z)b(z). 



(19) 



for a galaxy sample with the redshift distribution p(z)dz. 

Fig. [T] shows the time evolution of the source fields, i.e. the 
growth function D+(a) for the density field, and the time deriva- 
tive of D+(a)/a for the gravitational potential. Additionally, the lin- 
earised bias evolution model Ah{a) = feii„(a) - bo is plotted, for a 
range of bias parameters bo- As pointed out by a number of authors, 
the bias decreases with time and tends towards unity, although in 
the linear bias model (for small galaxy overdensities), the bias can 
never reach a value of exactly unity today if the galaxies form an 
initially biased population. 



4 ANGULAR POWER SPECTRA 



3.2 Galaxy density 

The bias prescription b(x) relates the galaxy number density to the 
density field S, 



7 ■ 



r" dxp(z)^b(x)D^(x)S. 
Jo 



dx 



(17) 



p(z)dz is the redshift distribution of the surveyed galaxy sample, 
transformed in terms of the comoving distance x, and D+(x) is 
the growth function of the density field at the scale factor a(x). 
We us e the redshift distribution of the main galaxy sample of EU- 
CLID teefregier & the DUNE collaboratioij|2008() . which will ob- 
serve half of the sky out to redshifts of or der unity and w hich will 
be of particular use for iSW-obser vations llDouspis et al.l[20()^ . In 
the parameterisation proposed bv lSmail et al.l 1 1995I) . the redshift 
distribution is approximated by 



p(z)d: 



■' = Po — I exp 



Zo 



with l = ^rp 

Po P XP 



(18) 



In summary, the line of sight integrals for the iSW-temperature per- 
turbation T and the galaxy density y read: 



3a„ r" , 2„, , d ID,\ 
nXH 

dxp{z)—b(x)DAx)S, 
Jo <iX 



(20) 
(21) 



where we have defined the dimensionless potential ip = A"'5/rf^, 
rescaled with the square of the Hubble distance = c/Ho for 
convenience. The weighting functions 



Wyix) 



a H{a)- , 

c da a 

p(z)^b(x)D4x), 



(22) 
(23) 



can be identified, which allow the expressions for the angular cross 
spectra to be written in a compact notation, applying a Limber- 
projection ( lLimbed[l954t) in the flat-sky approximation, for sim- 
plicity: 
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30 100 
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Figure 2. Angular iSW-spectra Cry{() (thick lines) and Cry{i)lb (thin lines) 
for the EUCLID main galaxy sample, for bias evolution parameters ba = 
(corresponding to the average bias b = \, solid line), ba = I (h = 1.46, 
dashed line), h„ = 2 (b = 1.93, dash-dotted line) and ba = 3 (h = 2.39, 
dotted line), for the EUCLID main galaxy sample with unit bias bg = i 
today. 



Jo 



iX — Ps^k 



with the cross-spectrum Ps^(k) 
spectra are given by: 



C/X) (24) 
Pd6(k)/(dHk)^. The angular auto- 



Jo 

/' 

Jo 



Ax 



dx 



W',(x) 

x' 



P^^ik = tlx). 



Pss(k = tlx). 



(25) 



(26) 



which will be needed in the expression for the covariance of Cry{t). 
In analogy to Ps^ik), the spectrum of the potential ip is defined as 
P^^ik) = P,,{k)KdHk)\ 

Angular cross power spectra Cry{() for a range of bias evo- 
lution paramters is depicted in Fig. |2] Assuming that an esti- 
mate of the average bias b follows from the spectrum Cyy{€) via 
b — CyyICry OX ftom other observational channels, the spectrum 
is rescaled with b because to first order, the iSW-spectrum Cry{{) 
is proportional to the mean bias. Fig. |2] suggests that on large an- 
gular scales, the difference in amplitude amounts to 7 . . . 12% for 
the range of b„ values considered here, relative to the constant bias 
model. This difference decreases towards smaller angular scales, 
which can be easily understood by decomposing H'^Ot') according 
to 



Wr(x)Wy(x) = W'(x)b(x) = W'Cr) [*o + (1 - a)ba] . 
The iSW-spectrum then consists of two contributions. 



(27) 



Cryit) = b 



Jo 



•XH W'ix) 

dx^^Ps, + b. 
X 



dx^Y^{l-a)P,^,(2%) 

Jo 



X 



(29) 



in analogy to the mean galaxy bias b, 

b = j dz p(z)b(z) = bo + b„ j dzp(z)(l-a), 

which is scale-independent and suppresses the amplitude in the 
scaled spectrum Cry{()lb, the suppression being proportional to 
balbo- The scaled angular spectrum, however, has in addition a 
scale-dependent contribution of the form 



t>o Jo 



X^ 



■ a)Ps^(k = eix). 



(30) 



whose integrand is large at early times, i.e. small a, which corre- 
sponds to large comoving distances x- Consequently, CDM spec- 
trum amplitudes corresponding to larger wave vectors k get pro- 
jected out, which leads to an overall decrease in the suppression, 
because Ps^ik) oc P(k)lk' is a monotonically decreasing function 
for the choice of the CDM transfer function T(k). The unsealed 
spectra Cryit), on the contrary, vary by almost half an order of mag- 
nitude, indicating the relative importance of bo and ba- 



5 FISHER-MATRIX ANALYSIS 

5.1 Fisher-matrix for the iSW-spectrum Cry(t) 

The Fisher matrix, which quantifies the decrease in likelihood if 
a model parameter moves away from the fiducial value, can be 
computed for a local Gaussian approximation to likelihood X 
exp(-^^/2). The Fisher-matrix for the measurement of Cjyit) is 
given by 



dCry(t) ,.^^^C,,(^) 
COV '\Cry{e),Cry{t)) 



dx. 



dxy 



(31) 



We construct the Fisher-matrix F^,,, for ACDM as the fiducial cos- 
mological model, where the parameter space is spanned by CI,,,, 
(Tg, h, n„ w and b„. The fiducial values for the parameters are 
fi„ = 0.25, CTg = 0.8, h = 0.72, n, = l,w = -1 and ba = 0, i.e. 
a non-evolving bias. It should be emphasised at this point that in a 
measurement of Cryit!) alone with ba = 0, the parameters as and bo 
are completely degenerate (leading to a singular Fisher-matrix), 

dCryjt) 2bo dCry({) ^^2^ 

das as dbo 
which is the reason why we set bo = 1, do not consider it as an 
independent parameter in the Fisher-analysis and absorb it into the 
normalisation of the spectrum as. 

Implicitly, we assume priors on spatial flatness, Q„, + fl^ = 1 
and on the primoridal slope of the CDM spectrum «, = 1, and ad- 
ditionally neglect the weak dependence of the shape parameter on 
the baryon density £1^. We extend the computation from f,^i„ = 10 
to ^max = 200, because at { t {^]„ the small-angle approxima- 
tion ceases to be applicable and because iSW-contributions from 
t ~ tmax are small and would be dominated by the nonlinear iSW- 
effect. Fig.[3]illustrates the sensitivity of the iSW-spectrum Cjyit) 
on variations of the cosmological parameters, by plotting the log- 
arithmic derivatives dlnCry{f)ld\nx^ (only in the case of the bias 
evolution parameter ba we give d\nCryi{)ldba instead as the ref- 
erence value ba = 0). One immediately recognises the degeneracy 
between erg and ba, which increases the mean bias, as well as be- 
tween f2,„ and h as they determine the CDM shape parameter. 

5.2 Noise sources and covariance 

In an actual observation, the iSW-power spectrum is modified by 
the intrinsic CMB-fluctuations, the instrumental noise and the beam 
as a noise source, assuming mutual uncorrelatedness of the individ- 
ual noise sources. The galaxy correlation function assumes a Pois- 
sonian noise term. 



\At) 

'yy(e) 



CrAt) + CcMB(i)+W-r' 

Cyy(£) + -, 



B-\£), 



(33) 
(34) 
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For PLANCK'S noise levels the value w^' = (0.02[iK)^ has 
been used, and the beam was assumed to be Gaussian, B~^(£) — 
exp(A02 ({t + 1)), with a FWHM- width of Ae = 7! 1, corresponding 
to the V = 143 GHz channel closest to the CMB-maximum. n in 
eqn. J34b corresponds to the number of galaxies per steradian. EU- 
CLID is expected to survey the entire extragalactic sky and to cover 
the solid angle Afl = 2n, i.e. /si^y = 0.5, yielding 4.7 x 10^ galax- 
ies per steradian. The observed cross power spectra are unbiased 
estimates of the actual spectra. 



(35) 



in the case of uncorrelated noise terms. For the determination of the 
spectrum Ccmb (^) of the in trinsic CMB anisotropies, the CAMB 
code written bv lLewis et al] JioOQ) was employed. The covariance 
of the spectrum Cry{{) is given in terms of the observed spectra 
Cxrit), with X,Y e jr, y), which follow directly from applying the 
Wick-theorem, 



Cov(C„,C„) : 



(36) 



with a cosmic variance error cx l//sky. Due to the relatively 
small signal to noise ratio attainable by iSW-measurements (the 
largest contribution to the covariance being the primary CMB 
anisotropies) we enhance the measurement by CMB priors. 



riSW 



(37) 



because Cry{t) and Ccmb(0 constitute independent measurements. 
The Fisher matrix describing the CMB parameter bounds is 
marginalised over the optical depth t and the baryon density 



5.3 Parameter bounds and degeneracies 

The ^--function for a pair of parameters {x^,Xy) can be computed 
from the inverse (-F"')^,, of the Fisher matrix. 



Ax^ 
Ax,, 



where Ax„ 



Ax,, 
Ax„ 



(38) 



two parameters is defined as 



The correlation coefficient r,,y between 



^{F-\,{F-\. 



(39) 



and describes the degree of dependence between the parameters x,, 
and Xy by assuming numerical values close to for independent, 
and close to unity for strongly dependent parameters. 

The degeneracies between the cosmological parameters Q„,, 
CTg, h, n,, w the bias evolution parameter is shown in Fig. |4] 
Most importantly, there is an anticorrelation between and erg, 
which both increase the iSW-signal Cry(l), likewise for ii, and w. 

Figs. [5] and m summarise the degeneracy of the cosmological 
parameters f2,„ and erg (i.e. the strength of the gravitational poten- 
tials) and of f2,„ and w (describing the particular dark energy model) 
on the bias evolution parameter ba- It is apparent in both cases that 
the inclusion of b^ as an additional parameter makes the parameter 
accuracy worse, and that the bias evolution b^ can be constrained 
to be < 0.35 by the iSW-data alone, with the CMB providing tight 
constraints on the cosmological parameters. 

The errors on the individual parameters, especially the 
bias evolution parameter ba, from a combined iSW- and CMB- 
measurement are given by Table [T] These errors follow from the 
Cramer-Rao bound: 




multipole order £ 



Figure 3. Logarithmic derivative d\nCry(t)ldx^ of the iSW-spectrum with 
respect to the cosmological pai'ameters x,, , specifically n„ (thick solid fine), 
(Tg (thick dashed Une), h (tfiick dash-dotted fine), n, (thin sofid fine), w (thin 
dashed fine) and ba (tliin dash-dotted fine). 



Figure 5. Fisher matrix constraint from Cisw(^) and Ccmb (t) on the tripfet 
(£1,,,, CTg, foo). The isoprobabifity surfaces correspond to f cr and 2a. The vec- 
tor (50.m,Scri,6ba) indicates the parameter estimation bias. 



T, = Ax, = ^/(F% 



(40) 



6 PARAMETER ESTIMATION BIAS 

The question of parameter estimation bias in the presence 
of systematics was addressed in different contexts: reionisa- 
tion (jZ ahn et al. 2005), weak cosmic shear ( Huterer et al. 20(3^; 
lAmara & Refregieil l2007) and Sunyaev-Zeldovich contaminations 
to the CMB spectrum iTaburet et alj2009|). In this s ection, we adapt 
the formalism worked out by l lTaburet et al.ll2009l) , because it can 
be used in the regime of strong systematic contributions to the ob- 
servable, in order to quantify how the interpretation of the data in 
cosmological model without bias evolution impacts on the estima- 
tion of cosmological parameters, when in reality the bias is evolv- 
ing. We restrict the parameter space to the cosmological parameters 
Q.„„ CTg, /j, and w, while assuming the bias to evolve according 
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Figure 4. Constraints on the parameters D,,,, erg, /i, vt' and ba, from a combined iSW and CMB, using ACDM as the fiducial model. The ellipses con'espond 
to lo" . . . 3tr intervals. Additionally, the vectors ((5^, rjy) point from the fiducial cosmology to the new best-fit position and indicate the bias in the estimation of 
the cosmological parameters if tracer bias evolution is neglected and if one falsely assumes = instead of = 1 . The number in the upper right corner of 
each panel gives the coiTelation coefficient r,,i, to two digits. 




Figure 6. Fisher matrix constraint from CiswOT and CcmbCO on the triplet 
mt ^1 ^ct)- The isoprobability surfaces correspond to Icr and Icr . The vec- 
tor i60.fn,Sw,^ba) indicates the parameter estimation bias. 



fiducial model 


statistical en'or 


estimation bias 


q 


n„ = 0.25 


An,,, = 0.038 


5n,„ = 0.045 


q=1.19 


o-g = 0.8 


Ao-g = 0.007 


<5o-g = -0.002 


q=0.28 


h = 0.72 


A/i = 0.021 


Sh = -0.022 


q=1.07 


"s = 1.0 


A(ij = 0.008 


Sn, = 0.005 


q=0.65 


w = -1.0 


Aw = 0.072 


Sw = 0.052 


q=0.72 


ba = 0.0 


Sba = 0.824 


Sba = 1.0 


q=1.21 



Table 1. Cramer-Rao bounds A.v^ on the parameter set as well as param- 
eter estimation biases 6x^ and the fraction g = |(5,,/Ap| between systematical 
and statistical error, from a combined CMB- and iSW-measurement, with 
ACDM and constant bias as the cosmological model. 



to b{a) = 2 - a (i.e. fco = 1 and ba = I), and while interpreting the 
data with a constant bias model b{a) = bg, b^ = 0. 

With our choice of the particular bias evolution model b(a) = 
bg + {I - a)ba we can decompose the distance weighting functions 



(41) 



such that the iSW cross-spectrum Cry{{) = Cisw(f) + Cevo(^) can be 
separated into a part with non-evolving bias Cisw(^) and an additive 
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systematic, containing tlie bias evolution C„oiO, by identifying tlie 
two contributions 



Jo 



(42) 



X 

C,M = K r" Ax^^[l-a{x)\P,^{k = llx). (43) 
Jo X 

Using these relations, we define the power spectrum of the true 
model C[{t) with evolving bias. 



C,(^) = CiSw(0 + Ce™(^) 



(44) 



as well as the spectrum of the false model C2{{), which neglects 
bias evolution. 



C2V) = 



(45) 



where the observed spectra Cj(i) are unbiased estimators of the the- 
oretical spectra Ci(£) in each case, because of the cross-correlation 
measurement method. 

If the data in reality is described by Ci{£), but is fitted erro- 
neously with Ciit), the best-fit values are biased, because the re- 
maining parameters have to emulate the ignored degree of freedom, 
in our case the bias evolution b^. This parameter estimation bias, 
defined as the distance 5 = X2-x\ between the best-fit values xi of 
the true model Ci(i) a nd X2 of the fal se model Ciit) can be derived 
using the formalism bv lTaburet et alj (2009), who expanded the;^'^- 
function of C2{t) at Xy in a Taylor series and recovered the best-fit 
position X2 by extremisation of the ensemble-averaged (x^): 



\dx. 



dXfjdx, 



■X 



<5„. 



(46) 



The resulting linear system of equations YivGi^ySy = can be in- 
verted for the estimation bias S, = 2v(G"');/rav The two quan- 
tities G^y and follow from the derivatives of the -function of 
model €2(0, evaluated atxi, 



fmax 



2Cov- 



fmax 



2Cov- 



flx^ dxy 
3C,sw(^) 



dxijdxy 



dxu 



(47) 



As in the case of the Fisher-analysis, the iSW-measurement is en- 
hanced by using CMB priors on the relevant parameters. A short 
calculation shows that the prior information can be in corporated 
by the replacement 

G,„ = G;^7 + ^^r, (48) 

in the case of uncorrelated iSW- and CMB-likelihoods. 

The biases in parameter estimation from the iSW-effect are 
depicted alongside the degeneracies in Fig. |4l for a bias evolution 
model with = I - Table[T]lists the parameter biases from the com- 
bined CMB and iSW-measurement: The shifts of the best-fit posi- 
tion are most prominent for the parameters as, iig and w, and are 
small for Q.,„ and h. Generally, the estimation bias is roughly oc b„ 
and the inclusion of a strong prior, like in our case of -F™^ , reduces 
the estimation bias significantly. The parameter bias is comparable 
to the statistical uncertainty at the Icr-level for the parameters erg, 
«, and w, which indicates that tracer bias evolution is an important 
systematic. Another point worth noting is that the estimation bias 
is not necessarily linked to the degeneracy direction. 

The estimation bias on w is particularly troublesome as w 
does not come out negative enough, and one would conclude from 



the measurement the presence of dark energy when in reality one 
would deal with a cosmological constant, erg is biased towards 
smaller values, because bias evolution implies a higher average 
bias in the past, such that erg can be lowered while retaining the 
same amplitude for the iSW-signal. Fitting the iSW-data with b^ 
as a free parameter deteriorates the statistical accuracy especially 
on Q.„, and w, which leads to the conclusion that it should be best 
to include external information on b^ from other data sets. Given 
the low significance of the iSW-signal, it is doubtful whether the 
data alone would prefer a more complicated model including bias 
evolution over a mod el with consta nt bias on grounds of Bayesian 
model selection CLiddle et alj2006l) . 



7 SUMMARY 

This paper treats the influence of bias evolution on the estima- 
tion of cosmological parameters from the iSW-effect, in particu- 
lar the dark energy equation of state parameter w. We consider 
a cross-correlation measurement with the galaxy density as bi- 
ased large-scale structure tracer and introduce two parameters bg, 
b„ as a generic bias evolution model, and assume that at red- 
shifts relevant for the iSW-effect, the bias is well described with 
b(a) = bo + (l-a)b,. 

(i) The iSW-spectra are sensitive to evolving bias model mainly 
through the amplitude of the mean bias of the tracer sample, with an 
additional scale-dependent contribution which is most appreciable 
on large angular scales. 

(ii) Parameter degeneracies between the cosmological parame- 
ters n„,, CTg, h, «, and w on one side and the bias evolution pa- 
rameter b„ on the other side were investigated in a classical Fisher- 
matrix approach, with ACDM as the fiducial cosmology. Due to 
the weakness of the iSW-effect, we consider combined constraints 
from the iSW-effect and from primary CMB fluctuations. We would 
like to emphasise that it is not possible to separate bias evolution 
from parameter estimation, but on the other hand find that b„ and w 
can be constraint simultaneously from the iSW-spectrum. 

(iii) In quantifying biases in the estimation of cosmological pa- 
rameters from the iSW-cross spectrum Cryit) if the data is inter- 
preted in a model with constant galaxy bias where in reality the 
bias is evolving we find that bias evolution shifts the best fit po- 
sition by an amount comparable to the statistical accuracy in the 
case of CTg, and w, indicating the importance of bias evolution on 
iSW parameter estimation. Specifically, the ratio between system- 
atical and statistical error amounts to 1.19 in f2„,, 0.28 in erg, and to 
0.72 in w for a bias evolution parameter of fc„ = 1. 

(iv) The parameters most affected are erg, due to the higher av- 
erage bias at earlier times, n, due to the change in the shape of the 
cross-correlation spectrum and w, which is not as well constrained 
as e.g. f2„, from primary CMB anisotropies. Ignoring bias evolution 
yields less negative values for the equation of state w, hinting at the 
danger that one might favour dark energy models over ACDM. 

In a future paper we will investigate means to control bias evolution 
by carrying out iSW-tomography as well as by combining the iSW- 
signal with biased tracers such as galaxies with the weak lensing 
signal from the same galaxy sample as an unbiased tracer of the 
potential fluctuations. 
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